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Abstract 



Subjecting the SU{2) Yang-Mills system to azimuthal symmetries in both the x — y 
! and the z — t planes results in a residual subsystem described by a [/(I) Higgs like model 

with two complex scalar fields on the quarter plane. The resulting instantons are labeled 



Q ' by integers (m, ni,n2) with topological charges q = ^ [1 — (— l)™']nin2. Solutions are 

constructed numerically for m 
only the m = 1 instantons are 
instanton-antiinstanton lumps. 



, constructed numerically for m = 1,2,3 and a range of ni = n2 = n. It is found that 

' only the m = 1 instantons are self-dual, the m > 1 configurations describing composite 



Oh! 

^ 1 Introduction 



After the discovery of the BPST instantons by Belavin et al pLj, which are unit topological 



' charge spherically symmetric solutions to the SU{2) Yang-Mills (YM) equations, Witten |2] 
constructed higher charge axially symmetric instantons. Multi-instantons subject to no sym- 
metries were subsequently constructed by 't Hooft and by Jackiw et al jS] for the same system, 
while the most general instantons for arbitrary gauge group SU{N) were classified by Atiyah 
et al [4- All these instantons are solutions to the first order self-duality equations. 

Solutions to the second order Euler-Lagrange equations which are not self-dual, while inter- 
esting for their own sake, are of great physical importance. The most important such solutions 
are the zero topological charge instanton-antiinstantons, whose putative role has been stud- 
ied from the eariest stages of instanton physics, in particular in the construction of dilute (or 
otherwise) instanton gases jSl E] • To date however no such exact solutions were constructed, 
although the forces between an instanton and an antiinstanton were considered long ago for 
approximate field configurations. 

The existence of non self-dual instantons has been proved in |H], and for the SU{2) YM 
system in CHI- In the latter, the residual actions studied are one dimensional as a result 
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of the imposition of quadrupole symmetry in J^, and by applying "equivariant geometry" in 
|10j . However, the proof in |9j does not cover the existence of non self-dual instantons with 
topological charges ±1. 

In this Letter, we present both self-dual and non self-dual instantons of the SU{2) YM 
model, subjected to azimuthal symmetries both in the x — y and the z — t planes. The residual 
system in our case is an axially symmetric U{1) Higgs like model in the p — a quarter plane^, 
as distinct from Witten's case |2J where it is a U{1) Higgs model on the half plane {r,t). We 
find finite action solutions labeled by a triple of integers (m, ni,n2) with topological charges 
q = I [1 — (— l)™]nin2. We find that only for m = 1 are the instantons self-dual, all those 
labeled by m > 2 being non self-dual. All our solutions are constructed numerically. 



2 The model 



= TTT^Tr , with J^^, = d^A, - d,A^ + [A^, A,], (1) 



2.1 Imposition of symmetry and residual action 

The usual SU (2) Yang-Mills (YM) action density in 4 Euclidean dimensions 

1 

16^2 

will be subjected to two successive axial symmetries. The normalisation used in is chosen 
such that the spherically symmetric BPST instanton has unit topological charge, in the 
spherically symmetric limit of our Ansatz to be stated below. We denote the Euclidean four 
dimensional coordinates as = {x,y; z,t) = {xa',Xi), with a = 1,2 and z = 3,4, and use the 
following parametrisation 

Xa = r sin 9 Xa = pXa, Xi = r COS 9 Xi = axi , (2) 

where = = |xq,P + with the unit vectors appearing in Q parametrised as Xa = 
(cosyji, sini^i), Xj = {cos ip 2, sin ip 2), with < 6* < | spanning the quarter plane, and the two 
azimuthal angles < y^i < 27r and < < 27r. 

The SU{2) YM field JF^,^ will be subjected to two stages of symmetry, in the Xa and the 
Xi planes, in succession. The first stage of symmetry imposition is of cylindrical symmetry in 
the Xa = {xi,X2) plane. Our cylindrically symmetric Ansatz is 

A = (^^y^) ^cpip + (^) (^^)" (^^^'^) + /3 - Al Xa4^ S^3 (3) 



A, = -4i4'^S^3 + ^-4'^S^4-4'S34 , (4) 

in terms of the unit vector ri'a' = (cosniy^i , sinni</)i) labeled by the vorticity integer ni, Sap 
being the Levi-Civita symbol. The spin matrices T^^y = (Sq,^, E^j, S^) in (jH)) and (jH) are one 
or other of the two chiral representations of 5*0(4), i.e. they are SU{2) matrices. 



^With = + and cr^ = + t^, where a;, z and t are rectangular cordinates on E' 
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Using the notation 



Fmn = dMA% - d^Alj + e'''"'AijA% , = du^ + e'^^'Al,^'^ , (5) 

with xm = {xi, p) and the internal index a = 1,2,3, the components of the curvature J-'^i, 
is expressed exclusively in terms of the gauge covariant quantities (jSl), describing an effective 
SO (3) YM-Higgs system on the hyperbolic space with coordinates xm = (yXi,p)- 

Imposition of the second stage of symmetry is in the Xj-plane in this hyperbolic space, by 
subjecting the fields {A'},A'^) and 0" appearing in (jS)), to axial symmetry in the Xi = (x^jX^) 
plane. Relabeling the S0{3) internal index as a = («',3), i.e. .. = 1,2, the axially 

symmetric Ansatz for the residual 5*0(3) gauge connection is 

A'- = -a^Xi (£72(2))*' + 1^^^ 71(2) , 

A^ - (^) , (6) 

A'p = -ap (671(2))'' , 

4 = 0, 

with the unit vector n'^^^^) = (cosn2V52, sinn2V92) labeled by a second vorticity integer n2. Sub- 
jecting the 5*0(3) triplet effective Higgs field 0" = (0*',0'^) to the same symmetry, we have 

In the numerical computation, we implemented a third stage of symmetry by treating the 
two azimuthal symmetries imposed in the x — y and the z — t planes on the same footing. 
Thus we set the two vorticities rii and n2, appearing in ©-© and (jH)) respectively, equal, 
rii = n2 = n. For the rest of this section however we consider the general case with the two 
distinct vorticities. 

As a result of imposing the two stages of symmetry, the YM field is parametrized by six 
functions which depend only on p and a. The components of the curvature (Fj",Fjp and the 
covariant derivatives {DpCp"-, Di(f)°-) appearing in (0) are now expressed exclusively in terms of 
the 5*0(2) curvature 

fpa = dptta - d^ap (7) 

and the covariant derivatives 

V,x'' = dpx'' + ap{ex)^ , V^x"" = d.x^ + a^{e'')^ . (8) 

Vpe = dpe+cip{eo'' , v^e = d.e+ciM)^. 

where we have used the notation (x^, x^) = X^^ {i^i ^^) = ^"^j A= 1,2. The result is a residual 
f/(l) connection (op, aa) interacting covariantly with the scalar fields x"^ and i.e. an Abelian 
Higgs like model in the quarter plane p — cr, described by a Lagrangean 



per fpa + ^ iiVpx^i' + iv.x^n + - imr + m^n + -(^^ve^)^ 

cr p 'pa 



(9) 
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This residual action density is a scalar with respect to the local SO (2) indices A, B, hence they 
are manifestly gauge invariant. It describes a U{1) Higgs like model with two effective Higgs 
fields ^-iid coupled minimally to the U{1) gauge connection = {ap,a^). To remove 
this U{1) gauge freedom we impose the usual gauge condition 

df, = dpttp + d„a„ = . (10) 

To state compactly the residual self-duality equations which saturate the topological lower 
bound of 0, we use a new index notation = (p, cr), not to be confused with the index 
notation used in (0). They are expressed in this notation as 

= -e,.e^\''e . (11) 
per 

V^X'' = p^'ae,,V,C- (12) 

One solution of these equations is well known, corresponding to the spherically symmetric 
unit charge BPST P instanton. Expressing the self-duality equations (fTT|l - (fT^ in terms of the 
coordinates (r, 6) rather than (p, a), constraining the function ag = ag{r) to be a radial function 
and assigning the following values of the remaining fuctions {arjX"^^^^)'- 

ar = 0, x^ = = sin 26*, x^ = ~'^e cos^ — 1, i'^ = — cte sin^ 6 — 1 , (13) 
these reduce to the single self-duality equation dag/dr = ag{ag + 2)/r yielding 

(A being the arbitrary scale of the unit charge instanton). 

Since our numerical constructions will be carried out using the coordinates (r, 6) we display 
(EI) also as 



rsin6' /.^ .,9 1 

+ - 

cos U \ 

rcos9 /,„ ^.,0 1 ,„ ^A,')\ 1 



rsin^cos^/2, + (\VrX^\' + -,\VgX^\' ] (15) 



1 = 1 

sm 6 \ J r sm U cos 6 

the total action of these solutions being S = f d^x^JgC = f^^'^ drdOL, with C the action 
density (jT)). 

It is worth noting that precisely the same results for the residual gauge connection and 
reduced action are found by the alternative approach of ^T] , where the action of the isometrics 
d/dipi, d/dip2 on the gauge connection is compensated by suitable gauge transformations. 
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2.2 Boundary conditions 

To obtain regular solutions with finite action density we impose at the origin (r = 0) the 
boundary conditions 

a, = , a, = , = f M , = f M , (16) 



^2 / V 

which are requested by the analyticity of the ansatz. In order to find finite action solutions, we 
impose at infinity 

= , a. = -2™ , = (-D-'+'n, ( ^^^^^ j , = -n. ( ^^^^^ j , (17) 

m being a positive integer. Similar considerations lead to the following boundary conditions 
on the p and a axes: 

a, = —dre , ae = —dee , = , = 0, = , = -r^i , (18) 

for = and 

a. = — , ae = —dex^ , = , = 0, = -^2 , 9^^' = , (19) 

^2 n2 

for Q = 7t/2, respectively. 

2.3 Topological charge 

Having determined the values of all the functions on the boundaries of (r, 6')-domain, we proceed 
to calculate the resulting topological charges. In our normalisation, the topological charge is 
defined as ^ 

1 or> ^^uupa I Tr^ au-F OCT d X ^ (20) 
OZTT^ J 

which after integration of the azimuthal angles (y^i, (^2) reduces to 

q = \j [sABX^'i'' fp. + Vypx^ V^^e) dp da , (21) 
which in the index notation used in (fTT|) - (fT^ can be expressed compactly as 

q = \e,u j {^-EABX^'e fpu + V.x^V^i''^ d^x (22) 

= \j e,. d,{x^V,i^ - e^uX^) d\ . (23) 

The integration in ()22|1 is carried out over the 2 dimensional ). As expected 

this is a total divergence expressed by (^Hj) . 
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Figure 1. The action density L is shown for the m 
and a = r cos 0\ 



1, n = 4 self-dual instanton (with p = rsin( 



Using Stokes' theorem, the two dimensional integral of (j23|) reduces to the one dimensional 
line integral 

q = ] [ X^V^^^ds^, (24) 



the line integral being taken around the loop 




V. 



p=0 



Ve e 



de + 



V, e 



dp . (25) 



The integral over the large quarter circle vanishes since Vqx^ and T>q^^ both vanish as r — » oo, 
so the the only contributions come from the a and p integrations. These are immediately 
evaluated by reading off the appropriate values of x^ ^"^^ from |T8|l - |T9| . The result is 



'-V. 



\n1n2 



(26) 



such that only for odd m is the Pontryagin charge nonzero and is then always equal to nin2- 
This is consistent with the description of instanton-antiinstanton chains. Obviously, if zero 
topological charge solutions exist, they must be non self-dual, as also must the higher odd 
m > 3 solutions solutions with nonzero charge since their actions are likely to grow with m 
while the (nontrivial) lower bound stays the same at nin2- It is likely that the odd m = 1 
solution does saturate this bound and hence is self-dual. We have verified all these features in 
our numerical constructions below. 



3 Numerical results 

Subject to the above boundary conditions (fTH|) - (fTm) we solve numerically the set of six coupled 
non-linear elliptic partial differential equations. We employed a compactified radial coordinate 
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Figure 2. The effective Lagrangean (a) and the topological charge density (b) as given by H15|l 
and (|22() respectively, are plotted for the m = 2, n = 2 instanton-antiinstanton solution. 

X = r/{l + r) in our computations, the equations being discretised on a nonequidistant grid in 
r and 6 with typical grid size 130 x 60. The numerical calculations were performed with the 
software package FIDISOL, based on the Newton- Raphson method [T^ . 

To simplify the general picture we set rii = n2 = n. We have studied solutions for m = 1 
with 1 < n < 5 and for m = 2 with 1 < n < 3. Our preliminary numerical results indicate that 
there exist also solutions with m = 3, and hence likely for all (m, n). 

Since the m = 1 solutions are special, we start by discussing their properties. It turns out 
that all solutions with m = 1 are self-dual. This was verified by checking numerically that the 
solutions of the second order equations satisfy the first order self-duality equations (fTT |) -|T ^ . 



7 



L- 
-2 



1 

logio(i') 



2 



3 



2 



1.5 



0.5 















\ / e=o \ : 








\ i e=m \ : 




\ i e=Tz/2 \ ; 



-2-1 1 2 3 

logio(r) 



Figure 3. The moduli |x| and |^| of the effective Higgs fields of the m = 2, n = 2 solution are 
shown as a function of r for several angles. 

Moreover the [m = l,n = 1) solution is spherically symmetric, corresponding to the unit 
charge BPST instanton (fT^ . (fT^ . The m = 1 solutions with n> 2 are axially symmetric and 
their gauge potentials a^, ae, X^i have nontrivial ^-dependence. A three dimensional plot 
of the action density for the m = 1, n = 4 self-dual instanton is presented in Figure 1. 

The m > 1 configurations satisfy only the second order Euler-lagrange field equations and 
are not self-dual. The function ag does not exhibit a strong angular dependence, while xi 
and ^1 have rather similar shapes. The package FIDISOL provides an error estimate for each 
unknown function. The typical numerical error is estimated to be on the order of 10~^, except 
for the n = 1 solutions which are somehow special. 
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Figure 4. Same as Figure 2 for the m = 3, n = 2 solution. 

In this case, although the numerical iteration still converges, the error is larger. This seems 
to originate in the behaviour of the function a^. which has a rather small ^—dependence (as 
opposed to the n > 1 configurations) and takes very small values. Its maximal error is around 
4% and comes from the vicinity of the origin. Thus our numerical results in this case are less 
conclusive, the existence and the properties of the m > 1 solutions with n — 1 requiring further 
work. 

A general feature of the solutions we found is that the action density C (or the Lagrangean 
L) posseses m maxima on the 6* = 7r/4 axis. Thus, in all cases studied it is possible to distinguish 
m individual concentrations of action, the relative distance between them being fixed by the 
location of the maximum values of the Lagrangean L (see Figures 2a and 4a) . 
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Figure 5. Same as Figure 3 for the m = 3, n = 2 solution. 

Increasing the value of n causes these maxima to become sharper but does not significantly 
affect the distances between the locations of these pseudoparticles. It is also interesting to note 
that, as shown in Figures 3 and 5 for the m = 2, 3 configurations with n = 2, the moduli of 
the effective Higgs fields |x| = ix'^X'^Y^'^ cind |^| = (^^^^)^'^^ possess m nodes on the p and a 
axis, respectively. The positions of these nodes coincides, within the numerical accuracy, with 
the locations of individual pseudoparticles. 

The topological charge density, namely the integrand of (|22j). also presents m local extrema 
on the 6 = tt/A axis, whose locations always coincide with the action density extrema. However, 
as shown in Figures 2b and 4b, the signs of the charges alternate between the locations of the 
successive lumps. 
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Thus the m = 1 configurations describe self-dual instanton solutions, m = 2 corresponds to 
an instanton-antiinstanton pair, while the m > 2 solutions are instanton-antiinstanton bound 
states composed of several pseudoparticles. 

The numerical results indicate that the action S{m, n) of a composite (m, n)— solution (with 
m > 1) is smaller than the action So{m,n) = mn^, of m single infinitely separated self- 
dual instantons each with vorticity n. We found e.g. S{2, 2) = 7.64, S{2, 3) = 15.96 while 
5(3,2) = 11.38. 

Also, we verified numerically by integrating ()22|) that, within the numerical accuracy, the 
solutions with m = 3 carry topological charge n^, while the topological charge of m = 2 
solutions vanishes for all values of n. 

4 Summary and discussion 

To summarise, we have constructed instantons of the four dimensional SU{2) YM system by 
numerically integrating the second order Euler-Lagrange equations of the residual two dimen- 
sional subsystem resulting from the imposition of azimuthal symmetries in both the x — y and 
the z — t planes. The residual system is a U{1) Higgs like model featuring two complex scalar 
functions. The instantons are labeled by a triple of integers [771,711, and have topological 
charges q = | [1 — {—l)"^]7iin2. Concrete constructions were presented for the cases m = 1, 2, 3 
and several values of 7ii = 712 = 71, although nontrivial solutions are likely to exist for any 
positive {771,711,712). 

Clearly, all solutions with even m carry vanishing topological charge, while all those with 
odd m carry charge 7ii 712 = 7i^, in our case. Thus our solutions describe instanton-antiinstanton 
lumps, rather analogous to the monopole-antimonopole chains of the Yang-Mills-Higgs model 

m 

With the exception of the instantons with m = 1, all the rest (with m > 2) do not saturate 
the topological bound and are only solutions to the second order field equations. They are 
non self-dual solutions with vorticity n, whose actions are always smaller than the action of m 
infinitely separated charge-n^, m = 1 self-dual instantons. This is supported by our numerical 
results. Instantons with m = 1 by contrast, do saturate the topological lower bound and satisfy 
the first order self-duality equations. Of these, the n = 1 solution is the charge- 1 spherically 
symmetric BPST instanton, while m = 1 instantons of charges are not spherically symmetric. 
These features are also borne out by our numerical results. 

The instantons and antiinstantons are located in alternating order on the p = a symmetry 
axis, at the nodes of the moduli Ix^^l and l^^^l of effective Higgs fields. A detailed study of these 
solutions will be presented in a future study. 

Concerning the relation of our results with those of CSl proving existence of non self-dual 
instantons, it may be interesting to speculate as follows: In the latter 01101 , it is stated that 
the proof of existence for non self-dual instantons of charge |g| = 1 is absent. In our case, while 
the self-dual solution with m = 1, tt, = 1 is evaluated in closed form, the m > 2, 71 = 1 non 
self-dual instantons of charge q = n"^ = 1 are the only solutions for which the numerical process 
had some difficulties, discussed in Section 3. This is the class of non self-dual instantons for 
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which the existence proofs of m > 2, n = 1 are absent. 

On a physical level, we hope that our solutions will prove useful both in the sense that they 
are labeled by integers (m, n) and describe higher charged instantons, and especially since they 
present exact non self-dual solutions that can be useful in the construction of instanton gases 
and liquids. Hopefully, these numerical results will be of help in constructing the solutions 
analytically. 
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